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I.  INTRODUCTION 


The  objective  of  the  present  contract  "as  to  develop  digital  tech¬ 
niques  for  the  processing  of  infrasonic  records.  The  records  consist 
of  infrasonic  signals  whose  source  is  remote  from  the  recording  stations. 

The  purpose  of  the  processing  is  to  obtain  the  best  possible  estimate  of 
the  infrasonic  signal. 

The  pure  infrasonic  signal  caused  by  an  event  is  modified  in  two 
different  ways  before  it  is  finally  recorded  at  locations  remote  from 
the  source:  (1)  as  the  signal  propagates  through  the  atmosphere,  meteoro¬ 
logical  parameters  such  as  wind  and  temperature  affect  the  signal,  and 
(2)  the  presence  of  noise  at  the  recording  sites  further  contaminates  the 
signal.  Thus  it  was  decided  to  develop  two  different  techniques  for  deal¬ 
ing  with  these  two  different  effects. 

To  eliminate  the  meteorological,  or  propagation,  effects  on  the  sig¬ 
nal,  an  empirical  study  was  begun  to  relate  the  observed  characteristics 
of  the  Infrasonic  signal  to  the  meteorological  conditions  prevailing  along 
the  path  of  propagation.  Unfortunately,  due  to  the  untimely  termination 
of  the  contract,  little  progress  was  made  in  this  area.  Such  a  study 
should  be  continued,  especially  in  view  of  the  fact  that  theoretical  work 
on  infrasonic  propagation  has  far  outpaced  sophisticated  observational 
studies . 

To  eliminate  or  suppress  the  noise  on  Infrasonic  records,  several 
digital  filtering  techniques  were  examined.  A  computer  program  that  uses 
a  prediction-error  filter  technique  was  developed  for  application  to  infra¬ 
sonic  records.  This  technique  uses  samples  of  noise  prior  to  the  infra¬ 
sonic  signal  to  develop  a  Wiener  prediction  filter,  which  when  applied  to 
the  infrasonic  record,  predicts  the  noise  during  the  start  of  the  signal. 
Subtraction  of  this  estimate  of  the  noise  from  the  original  record  provides 
an  improved  estimate  of  the  infrasonic  signal. 

A  short  discussion  of  our  planned  activities  on  the  empirical  studies 
of  atmospheric  propagation  effects  on  infrasonic  signals  is  contained  in 
Section  II  of  this  report.  Section  III  contains  a  oiscussion  of  the  several 
possible  digital  filtering  techniques  that  were  examined;  an  explanation 
of  the  prediction-error  filter  technique  that  was  adopted;  an  analysis  of 
the  results  obtained  from  application  of  the  prediction-error  filter  tech¬ 
nique  to  actual  infrasonic  seconds;  and  a  detailed  description  of  the  com¬ 
puter  program  that  was  developed  to  accomplish  the  filtering.  Section  IV 
outlines  our  conclusions  and  contains  suggestions  for  further  evaluation 
and  improvement  of  these  techniques. 
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II.  EMPIRICAL  STUDIES  OF  CONCURRENT  INFRASONIC 
AND  METEOROLOGICAL  DATA 


A  thorough  review  of  the  literature  on  the  Infrasonlcs  problem  was 
performed.  Most  of  the  work  lu  the  literature  concerns  the  development 
of  theoretical  models  to  predict  the  propagation  of  acoustic-gravity 
waves  in  the  atmosphere.  The  early  theoretical  models  dealt  solely  with 
the  propagation  problem  in  a  non- isothermal  atmosphere  and  did  not  in¬ 
clude  the  effect  of  winds  (see,  for  example,  Pfeffer,  1962).  Later 
models  included  modeling  of  the  explosive  source  (for  example,  Harkridor, 
1964),  but  still  did  not  include  winds.  The  more  recent  work  of  Mac¬ 
Kinnon  (1968)  and  Pierce  (1968)  includes  the  effect  of  winds  and  indi¬ 
cates  that  winds  are  extremely  important  in  the  determination  of  the 
characteristics  of  infrasonic  signals.  The  work  of  MacKinnon  shows  that 
the  winds  have  a  pronounced,  but  complicated  effect  upon  the  amplitudes 
of  the  infrasonic  signals.  Thus,  recent  theoretical  work  suggests  that 
in  order  to  arrive  at  estimates  of  source  characteristics  from  infrasonic 
records,  the  effect  of  the  atmosphere  on  the  signal  should  be  eliminated. 

In  our  review  of  the  literature,  we  found  few  observational  analyses 
that  attempted  r.o  relate  observed  characteristics  of  infrasonic  waves  to 
'.he  prevailing  meteorological  situation.  Furthermore,  the  few  that  were 
found  (Wexler  and  Hass,  1952;  MacKinnon,  1968)  used  rather  gross  analyses 
of  the  meteorological  situation. 

As  a  result  of  our  review  of  the  literature,  we  arrived  at  certain 
cone lusions. 

(1)  The  effect  of  the  atmosphere  on  the  characteristics  of  the  ob¬ 
served  infrasonic  signal  can  be  considerable. 

(2)  To  estimate  source  strength  from  infrasonic  signals  one  should 
eliminate  the  effect  of  the  atmosphere. 

(3)  Previous  empirical  analyses  of  infrasonic  data  and  concurrent 
meteorological  data  are  either  non-existent  or  were  grossly 
performed. 

Therefore,  we  planned  to  perform  empirical  analyses  of  infrasonic 
data  and  concurrent  meteorological  conditions.  Guided  by  the  results  of 
the  recent  theoretical  work,  one  should  be  able  to  develop  empirical  rela¬ 
tionships  between  the  meteorological  conditions  -  in  particular,  integrated 

great  circle  winds  and  temperatures  for  key  layers  in  the  atmosphere  -  and 
characteristics  of  the  signal,  such  as  amplitude  and  period.  The  goal 
would  be  the  development  of  an  objective  procedure  for  "normalizing"  for 
prevailing  atmospheric  conditions.  Given  such  a  procedure,  one  would  be 
in  a  much  better  position  to  estimate  source  characteristics  from  infra¬ 
sonic  signal  characteristics  measured  remotely  from  the  source. 
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For  this  purpose,  we  obtained, from  the  National  Weather  Records  Center 
at  Asheville,  microfilm  copies  of  analyzed  daily  surface  and  upper  air 
charts  for  the  Northern  Hemisphere  for  times  corresponding  to  infrasonic 
signals. 

Unfortunately,  due  to  the  untimely  termination  of  the  contract,  this 
research  task  was  not  carried  out.  The  meteorological  data  are  being  trans 
ferred  to  the  contract  monitor  under  separate  cover.  It  is  hoped  that  such 
a  study  can  be  conducted  in  the  near  future. 
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III.  DIGITAL  FILTER  STUDIES  TO  DEVELOP  A  NOISE 
AND  BACKGROUND  SUPPRESSION  SCHEME 


A.  Basic  Problem 


The  basic  problem  In  the  enhancement  of  Infrasonlc  records  can  be  sum> 
marized  as  follows.  Infrasonlc  records  are  available  from  several  instru¬ 
ments  (channels)  at  slightly  different  locations.  A  schematic  diagram  of 
such  records  from  a  3  channel  network  is  shown  as  Figure  1.  The  record 
consists  of  a  period  of  noise  followed  by  a  period  of  time  in  which  the 
record  consists  of  noise  plus  an  infrasonlc  signal.  Aside  from  possible 
time  lags,  the  signal  should  be  exactly  the  same  at  all  three  channels. 

The  major  characteristic  of  the  signal  is  an  initial  high  amplitude  oscil¬ 
lation,  which  we  shall  call  the  first  motion.  This  is  generally  followed 
by  smaller  amplitude  oscillations.  The  problem  is  to  take  the  records  from 
the  different  channels  and  digitially  manipulate  them  to  obtain  the  best 
estimate  of  the  infrasonlc  signal. 

B.  Digital  Filters 

To  solve  the  problem  of  obtaining  the  best  estimate  of  an  infrasonlc 
signal,  we  have  reviewed  several  different  types  of  di,5ltal  filtering 
techniques : 

(1)  Wiener  filter 

(2)  Maximum  likelihood  filter 

(3)  Prediction-error  filter 

The  general  Wiener  filter  is  an  enhancement  and  prediction  filter. 
Enhancement  is  defined  as  separation  of  the  signal  from  the  signal  and 
noise  time  series.  Prediction  is  defined  as  the  forecast  of  the  signal 
for  times  greater  than  t,  based  upon  information  available  up  to  time  t. 

In  general,  one  has  a  signal  and  noise  time  series,  >  which  consists  of 
a  signal  s^,  and  noise,  n  ,  so  that 


=  s^  + 


(1) 


The  goal  is  to  design  a  physically  realizable  (physically  realizable  means 
that  the  filter  uses  only  the  past  history  of  the  time  series)  digital  filter 
so  that  when  it  is  applied  to  xj.,  the  actual  output  y^  is  the  best  approxi¬ 
mation  to  the  desired  output  Zj.  =  sj.  +  q  where  a  represents  the  prediction 
span.  In  roost  cases  of  Wiener  filtering,  it  is  assumed  that  the  form  of 
the  signal  is  known,  and  that  this  is  the  desired  output.  The  criterion 
employed  in  the  derivation  of  the  filter  is  the  minimization  of  the  squared 
difference  between  the  filtered  time  series  (output)  y  and  the  estimate 


CHANNEL  1 


IRST  MOTION 


Schematic  diagram  of  infrasonic 
records. 


of  what  the  time  series  would  be  if  noise  were  absent  (desired  output)  z^. 
Unfortunately,  in  the  case  of  the  infrasonic  records  we  do  not  know  pre¬ 
cisely  the  form  of  the  signal  and,  hence,  we  do  not  use  Wiener  filtering 
directly. 

The  maximum  likelihood  filter  (see,  for  example,  Green,  et  al.,  1966) 
maximizes  output  signal-to-noise  ratio  under  the  restriction  of  zero  fre¬ 
quency  distortion  of  the  signal.  It  has  been  utilized  in  the  analysis 
of  seismic  data,  where  it  is  derived  in  the  following  manner.  A  sample 
of  the  noise  record  preceding  an  event  is  used  to  derive  filter  coefficients 
which, when  applied  to  the  noise  record,will  produce  an  output  record  of 
minimum  variance  (that  is,  will  tend  to  drive  the  record  to  zeros).  When 
this  filter  is  passed  over  the  entire  length  of  record  (that  is,  the  noise 
portion  plus  the  portion  containing  signal  and  noise),  the  record  is  driven 
to  zero  in  the  noise  portion  and  to  the  signal  in  the  signal  and  noise 
portion.  However,  this  is  not  a  predictive  filter;  it  is  implicitly  assumed 
that  the  characteristics  of  the  noise  do  not  change  in  the  signal  and  noise 
portion  of  the  record.  It  is  of  possible  application  to  the  infrasonics 
problem  since  the  only  requirement  for  its  development  is  the  availability 
of  a  sample  of  noise  preceding  an  event.  As  opposed  to  the  Wiener  tech¬ 
nique,  one  need  know  nothing  about  the  characteristics  of  the  signal. 

The  prediction-error  filter  (see,  for  example,  Claerbout,  1964)  com¬ 
bines  features  of  both  the  Wiener  filter  and  the  maximum  likelihood  filter. 

A  prediction-error  filter  predicts  the  noise  at  a  future  time.  It  thus 
also  differs  from  both  the  Wiener  filter,  which  attempts  to  predict  the 
signal  at  a  future  time,  and  the  maximum  likelihood  filter,  which  attempts 
to  estimate  the  noise  at  the  present  time.  With  a  prediction  of  the  noise, 
it  is  possible  to  subtract  the  predicted  noise  from  the  actual  record,  thus 
obtaining  an  improved  estimate  of  the  signal.  The  prediction-error  filter 
will  obviously  work  well  if  the  noise  has  some  degree  of  predictability. 

Like  the  maximum  likelihood  technique,  it  is  applicable  to  the  infrasonics 
problem,  and, in  principle,  both  should  yield  similar  results.  Because  of 
the  availability  of  a  number  of  computer  sub-routines  (see  Robinson,  1967) 
that  are  directly  appllcable,we  decided  to  test  this  filtering  technique 
on  the  infrasonic  records,  ^e  technique  is  described  in  the  next  section 
(C);  the  computer  program  is  presented  in  Section  D;  and  the  results  are 
discussed  in  Section  E.  Further  discussions  of  the  various  filtering  tech¬ 
niques  and  our  circuitous  route  to  development  of  the  final  filtering 
technique  are  contained  in  the  contract  quarterly  reports. 

C.  Prediction-Error  Filter  Technique  for  Analysis  of  Infrasonic  Records 

The  prediction-error  filter  forecasts  the  noise  at  a  future  tine.  In 
general,  the  ability  of  such  a  filter  to  predict  the  noise  will  decrease 
with  the  length  of  the  prediction  span  (time  Interval  between  time  of  last 
data  input  and  time  of  predicted  noise).  Ihus,  if  possible,  it  is  desirable 
to  have  the  prediction  span  as  short  as  possible  for  a  particular  problem. 

In  the  infrasonics  problem,  as  discussed  above,  a  major  feature  of  the 
signal  is  an  initial  relatively  large  amplitude  oscillation.  We  may  call 
this  oscillation  the  first  motion. 
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We  shall  assume  that  much  of  the  information  on  the  event  that  causes 
the  infrasonic  signal  is  concentrated  in  this  first  motion.  Therefore, 
if  we  can  obtain  a  good  estimate  of  the  first  motion,  we  shall  be  in 
good  shape.  This  is  obviously  easier  than  obtaining  a  good  estimate  of 
the  entire  length  of  signal,  which  requires  a  very  long  prediction  span. 
Thus,  for  the  analysis  of  Infrasonic  records  we  shall  use  a  prediction 
span  equal  to  the  time  period  during  which  the  first  motion  occurs. 

In  general,  a  prediction-error  filter  for  a  network  of  channels  will 
use  the  history  of  the  ith  channel  time  series  and  the  history  of  all  the 
other  time  series  in  the  network  to  predict  the  noise  on  the  ith  channel. 
When  this  noise  is  subtracted  from  the  1th  channel  during  the  noise  plus 
signal  portion  of  the  record,  an  estimate  of  the  signal  is  obtained. 


The  filter  is  derived  from  the  noise  samples  preceding  the  event,  and 
can  be  derived  as  a  Wiener  filter  that  uses  estimated  future  noise  rather 
than  estimated  future  signal  as  the  desired  output.  In  the  development 
below,  we  follow  Robinson  (1967).  Let 


=  x(t)  = 


X2(t) 


=  input  time  series  (samples  of 
noise  preceding  the  event) 


(2) 


=  2(t) 


X^(t)' 

X2(t) 


desired  output  time  series  (3) 

[  x(t  +  a)  ] 


Yj.  s  y(t)  = 


y^Ct) 

y2(t) 


=  output  time  series 


(4) 


where  N  is  the  number  of  input  channels,  M  is  the  number  of  output  chan¬ 
nels,  and  a  is  the  prediction  span.  The  output  time  series  is  obtained 
from  the  convolution 
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y  =  (f)(x) 


(5) 


where  £  is  Che  desired  filter,  f  Is  oerived  by  minimizing  the  expression 
-  y^)^.  The  desired  output  time  series  is  simply  the  noise  pre¬ 
ceding  Che  event  shifted  in  time  by  the  prediction  distance  a.  The  output 
time  series  y^  represents  predictions  of  the  noise  at  absolute  time  (t+a) 
based  upon  application  of  the  filter  £  to  the  input  series  at  time  t. 

The  Wiener  digital  filter  Co  predict  the  noise  is  based  upon  three 
assumptions  which  are  implicit  in  its  derivation  and  during  its  application: 

(1)  The  noise  statistics  remain  stationary  for  a  time  outside  the 
noise  sample  Interval. 

(2)  The  approximation  criterion  is  the  minimization  of  the  mean- 
square-error  matrix  between  the  desired  output  and  the  actual 
output. 

(3)  The  filters  are  assumed  to  be  linear  and  physically  realizable. 

The  success  of  the  filter  depends  essentially  upon  the  existence  of  correla¬ 
tions  between  the  desired  output  and  the  input  x^.  The  assumption  of 
linear  filters  implies  use  of  linear  correlation,  namely,  the  autocorrela¬ 
tion  function 

=  E  |xj.  x^_^|  ,  T  =  0,1,2 -  (6) 


of  Che  input,  and  the  cross-correlation  function 

(t)  =  E  fz^  x^  1  »  T  =  0,1,2.... 
zx  (,  t  t-x) 

between  desired  output  and  In^uc. 

The  filter  may  be  represented  by  the  M  x  N  matrix-valued  cofficients 


(7) 


f  ^  f(s)  = 

! 

s 

CD 

f •  • ' 

...  f,jj(s) 

(8) 


It  is  assumed  Chat  mean  values  have  been  removed  from  the  input  and  desired 
output,  so  E  {x^l  and  E  {z^}  =  0.  The  actual  output,  which  is  the  convolu¬ 
tion  of  the  input  with  the  filter,  is  given  by  the  matrix  equation 


f  x  +  f, 
o  t  1 


Vi-^ 


.+  f  X. 
m  t- 


m 


(9) 


where  m  is  the  length  of  the  filter. 
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The  error  e^  between  desired  output  and  actual  output  is  Chen 

®t  “  *t“^t  “  *t  "  ^^o  *t  *t-l  *t-m^ 

The  tnean-square-virror  matrix  is  defined  by  E  written 

as  ^  ^ 

E|ej^(t)e2(t)|....  E|ej^(t)  e^j(t)| 

^ ^  'K<'>  '{4  '4  _ 


(11) 


The  values  of  the  filter  coefficients  fg  are  determined  from  the  condition 
that  the  trace  (abbreviated  tr)  1  of  the  mean-square -error  matrix  be  a 
minimum.  This  is  accomplished  by  setting  Che  partial  derivative  of 


I  «  tr  E|Vt^}  '  E|e2^(t)J+ . +  E  (li) 


with  respect  to  each  filter  coefficient  equal  to  zero.  It  can  be  shown 
that  when  this  is  done,  the  following  normal  equations  are  obtained 


f  ♦  (0)  +  f,  o  (-1)  +  ...  +f  ♦  (-m)  »  ♦  (0) 

o  xx'  1  xx'  '  m  XX  '  '  zx  '  ' 

f  0  (1)  +  f,  O  (0)  +  ...  +f  ♦  (1-m)  -  ♦  (1) 

o  xx'  *  1  xx'  '  m  XX  '  '  zx' 


(13) 


f  4i  (m)  +  f,  <t>  (m-1)  +  ...+!♦  (0)  *  ^  (m) 

o  xx'  '  1  xx'  '  m  XX  '  '  zx' 


The  known  quantities  are  the  correlation  ccefficients  and  the  unknowns 
are  the  filters  (f  ,  f  . . . .  ,f^) ,  which  can  be  derived  by  solving  these 
equations.  Procedures ^for  solution  and  further  details  in  the  deriva¬ 
tion  are  given  by  Robinson  (1967). 


D.  Computer  Program  MAXLKH 

To  accomplish  the  desired  digital  filtering,  a  computer  program  en¬ 
titled  MAXUQl  was  developed  and  tested.  This  program  makes  use  of  a 
number  of  sub-routines  presented  in  Robinson  (1967).  Basically,  what 
MAXUCH  does  is  the  following.  It  reads  in  infrasonic  time  series  (x^) 
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for  a  number  of  channels.  The  first  portions  of  these  time  series  repre¬ 
sent  noise,  the  second  portions  represent  noise  plus  signal.  The  first 
portions  are  used  as  noise  samples,  are  displaced  to  the  left  in  time  by 
the  prediction  span  time,  and  become  the  desired  output  series  for 
Wiener  filter  derivations.  Wiener  filters  are  derived  and  are  applied  to 
the  infrasonic  records.  This  application  results  in  an  output  time  series 
for  each  channel  that  represents  estimates  of  the  noise  for  the  time  period 
that  covers  the  noise  portion  and  the  first  motion  portion  of  the  infra¬ 
sonic  record.  These  estimates  of  noise  are  subtracted  from  the  original 
infrasonic  record  to  obtain  a  predicted  time  series.  The  predicted  time 
series  should  have  minimal  noise  in  the  noise  portion  of  the  record  and  an 
enhanced  first  motion  in  the  noise  plus  signal  portion.  For  the  time  period 
after  the  first  motion,  the  predicted  time  series  is  assumed  to  be  the  same 
as  the  original  time  series.  Program  MAXLKH  includes  a  plotting  sub-routine 
which  permits  the  input,  assumed  output,  output,  and  predicted  time  series  to 
be  presented  on  graphs. 

The  remainder  of  this  section  contains  information  that  is  required  to 
understand  and  operate  this  computer  program.  A  flow  chart  of  Program 
MAXLKH  is  shown  in  Figure  2.  A  list  of  variables  appearing  in  MAXLKH  is 
contained  in  Table  1,  and  dimensions  of  arrays  contained  in  the  program 
are  shown  in  Table  2.  The  content  of  the  data  input  cards  is  shown  in 
Table  3. 

A  program  card  deck  is  being  supplied  to  the  contract  monitor  under 
separate  cover.  To  operate  the  program  one  must  check  those  cards  that 
hsve  a  blue  top.  These  cards  can  vary  from  case  to  case,  and  can  vary  if 
the  filter  length  is  changed. 


E.  Results  of  Application  of  Program  MAXLKH  to  Infrasonic  Records 

Traces  of  infrasonic  records  were  provided  by  the  contract  monitor. 
These  were  digitized  and  punched  on  cards  for  use  as  input  to  program 
MAXLKH.  Figure  3  Indicates  the  results  of  application  of  MAXLKH  to  infra¬ 
sonic  data.  Figure  3  is  divided  into  3  parts;  3a,  3b,  and  3c,  each  part 
representing  a  different  channel.  On  each  part  are  four  time  series.  From 
top  to  bottom,  they  are  the  input  time  series  (x^.,  the  original  infrasonic 
record),  the  assumed  output  time  series  (z^,  the  original  infrasonic  record 
shifted  to  the  left  by  an  amount  equal  to  the  prediction  span  «  ISKIP) ,  the 
actual  output  time  series  (y^,  predicted  noise  time  series  plus  zeros  before 
the  beginning  and  after  the  end  of  the  noise  series),  and  the  time  series 
representing  the  difference  between  the  input  time  series  and  the  output 
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Figure  2.  Flow  Chart  of  MAXLKH 


TABLE  1 


AVG(I) 

E(I) 

F(I,J,K) 

FLOOR 

FORMAT 

IC 

ID 

IP 

ISXIP 

L 

LENGTH 

LF 

LR 

LW 

LX 

LXl 

LY 

LZ 

LZSKIP 

LZl 

M 

N 

NLENG 


LIST  OF  VARIABLES  APPEARING  IN  PROGRAM  MAXLKH 

llie  average  value  of  the  I^^  channel. 

The  mean  square  error  for  a  filter  of  length  I. 

The  filter  for  an  I  channel  Input  by  J  channel  ovcput 
by  K  filter  length. 

Ignore 

Array  continuing  data  Input  format. 

Card  reader  device  number. 

Array  for  plotter  which  titles  first  frame  identification. 
Printer  device  code  number. 

Prediction  span. 

Length  of  predicted  noise  series  (Filter  length  + 
input  length  NLENG). 

Length  of  noise  sample  to  be  used  in  Wiener (=  NLENG  - 
ISKIP) . 

Length  of  filter 
Max.  length  of  filter 
Not  used 

Length  of  initial  input  series 

Length  of  noise  sample  used  i  finding  filters  =  LENGTH 
Length  of  output  (LX-ISKIP)  +  LF-1 
Length  of  assxmed  output  =*  LX-ISKIP 
LX  -  ISKIP 

Length  of  array  used  as  for  Wiener  program 
(length  -  LXl) 

Number  of  output  channels 
Number  of  input  channels 

Maximum  length  of  data  available  to  be  used  in  finding 
the  filters 
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TABLE  1  (Continued) 


S  (I) 
TITLE  (7) 

X  (I.J) 
XY  (I) 

XI  (l.J) 

Y  (I.J) 

Yi  (I.J) 

Z  (I.J) 

Z1  (I.J) 


Storage  array  for  auto  and  cross  correlations 
Array  which  contains  the  title  of  the  data  set 
Input  array 

Temporary  storage  array 

Array  of  data  used  to  find  the  filters  (length  =  LXl) 
Output  array 

Output  array  for  data  used  to  define  the  filters 
Assumed  output  Z(I)  =  X(I  +  ISKIP) 

Assumed  output  for  data  used  to  find  Wiener  filters 
(length  *  LZl) 
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TABLE  2 


AVG  (I) 

XY 

(I) 

Title  (7) 

Format  (10) 

ID 

(10) 

X 

(I.J) 

z 

(I,J) 

F 

(I,J,K) 

E 

(I) 

Y 

(I.J) 

S 

(I) 

XI 

(I,J) 

Z1 

(I.J) 

Y1 

(1,J) 

DIMENSIONS  OF  ARRAYS  APPEARING  IN  PROGRAM  MAXLKH 

Where  I  >  number  of  channel* 

Where  I  «  maximum  length  of  X  or  Z  series 

Fixed 

Fixed 

Fixed 

I  is  the  number  of  input  channels,  J  >  channel  length 
I  is  the  nimiber  of  output  channels,  J  «  channel  length 
I  *  number  of  output  channels,  J  =  number  of  input  channels, 
K  >>  maximum  filter  length 
I  maximum  filter  length 
I  ••  number  of  output  channels,  J  ••  LX  *f  LF>1 
I  =  (N*N*5*LR)  +  (M-N*LR)  +  (M*N) 

I  s  number  of  input  channels,  J  >  length  of  data  used 
I  »  number  of  output  channels,  J  »  LZl  LXl 
I  =  number  of  output  channels,  J  «  (LXl  +  LF-1) 


TABLE  3 


CONTENT  OF  DATA  INPUT  CARDS  FOR  PROGRAM  MAXLKH 


1®*"  card 
,2”**  card 

3*^*^  card 

4*"^  card 

„th  , 
N  card 


sti 

Title  of  data,  1  40  columns 

N,LX,M,LZ,LR,Uf,FLO(Xt,  ISKIP,  NLENG 
format  <6I5,F5.3,3I5) 

Data  format,  columns  1  to  80 
Input  data 
Final  input  data 
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Seq.  ^^5656AIL 

ISRIP  =  Prediction  Span  = 

150 

Filter  Length  = 

25 

NLENG  *  Noise  Saaple 

750 

\  ■ 

(ARRAY  X) 

ThESC  0**a  «ii|  iAOM  AuM  ttS4l 

THANNli  *<0  J 

—  ' 

LX-1800 • - 

OUTPUT  (ARRAY  Z) 


NLENG  -■ 


ISKIP-* 


LZ=LX-ISKIP=L2SKIP 


I  sirs 


9»« 

mm%  SmiW ttm  i 


Ml  ••  ■  Itl  •Mil  I 


loc  1^0 

LENGTH—.. . 

(ARRAY  Y) 


io«o  ‘Tiw 


Figure  3a.  Application  of  Program  MAXLKH  to  Infrasonic  Record  11343 
Channel  3. 


eq.  #5656AIL 


Figure  3c.  Channel  2 


ftrrtrrrni 


time  series  (x^.  -  y^).  This  last  time  series  represents  the  filtered 
time  series  and  should  have  minimum  noise  prior  to  the  event  and  enhanced 
first  motion  during  the  event. 

Figure  3k  also  shows  some  of  the  variables  used  In  the  program.  In 
this  case,  the  prediction  span  (ISKIP)  is  150  time  tmits  which  covers  the 
first  motion  and  some  additional  oscillations.  The  total  noise  sample 
prior  to  the  event  has  a  length  (NLENG)  equal  to  750  time  imits.  However, 
since  the  prediction  span  is  150  units,  not  all  of  the  total  noise  sample 
is  available  to  derive  the  Wiener  filter.  After  being  shifted  to  the  left 
150  unii.s,  the  length  of  noise  sample  actually  used  (LENGTH)  in  deriving 
the  filter  is  600  time  units.  These  var1.ous  lengths,  plus  several  others, 
are  illustrated  in  Figure  la. 

The  success  of  the  filter  can  be  evaluated  by  comparing  the  input  time 
series  (top)  with  the  filtered  time  series  (bottom).  For  all  three  channels  it 
can  be  seen  that  application  of  the  filter  has  resulted  in: 

(1)  A  decrease  in  the  noise  prior  to  the  event. 

(2)  An  increase  in  the  amplitude  of  the  first  motion. 

The  reduction  of  noise  prior  to  the  event  makes  the  first  motion  easier  to  iden¬ 
tify.  The  increase  in  a&q>litude  of  the  first  motion  is  presumably  due  to 
the  fact  that,  in  this  case,  noise  had  acted  to  reduce  the  amplitude  of  the 
first  motion.  The  resulting  first  motion  should  thus  represent  an  improved 
estimate  of  the  signal  prevailing  at  this  time. 

Application  of  the  filtering  technique  to  another  set  of  infrasonic 
records  is  illustrated  in  Figure  4.  In  this  case,  the  prediction  span  is 
somewhat  shorter,  100  time  units.  Comparison  of  the  filtered  time  series 
with  the  unflltered  series  reveals  a  large  reduction  in  noise  prior  to  the 
event  on  all  four  channels.  Presumably,  there  has  been  a  comparable  noise 
reduction  during  the  time  period  covered  by  the  first  motion.  Thus,  the 
filtered  time  series  should  contain  an  improved  estimate  of  the  first 
motion. 

If  the  filters  worked  perfectly,  the  first  motions  in  the  filtered 
time  series  would  be  identical  on  all  channels.  Since  they  are  not  per¬ 
fect,  the  first  motions  differ  somewhat.  An  improved  estimate  of  the 
true  value  of  the  first  motion  can  be  obtained  by  averaging  the  individual 
filtered  time  series  to  obtain  a  mean  value  of  the  filtered  time  series 
during  the  first  motion.  This  feature  is  not  presently  included  in  pro¬ 
gram  MAXLKH,  but  can  easily  be  added. 


A  last  example  of  the  application  of  program  MAXLKH  to  infrasonic 
data  is  shown  in  Figure  5.  In  this  case,  no  infrasonic  signal  is  present 
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Figure  4a.  Application  of  Program  MAXLKH  to  Infrasonic 
Record  11354;  Channel  1. 


See.  #5433ANI 


Prediction  Span  =  100 
Filter  Length  =  25 
Noise  Sample  =  300 


Prediction  Span 
Filter  Length  *= 
Noise  Sample  >» 


Seq.  #5433ANI 


100 

2'i 

300 


IHCSt  CaU  FitCM  lUK  .hkHALl  K  i 


Seq.  #5628ADV 


Prediction  Span  *  18 
Filter  Length  »=  18 
Noise  Sample  =  160 


Figure  5a.  Application  of  Program  MAXLKH  to  Infra.son.c 
Noise  from  Record  63010,  ('.liannel  1. 


Seq.  #5628ADV 


Prediction  Span  =  18 
Filter  Length  -  18 
Noise  Sample  3:  X60 


and  the  objective  of  the  experiment  is  to  determine  how  well  a  typical 
sample  of  noise  can  be  predicted.  The  prediction  span  is  short  >  13  time 
units.  A  Wiener  filter  of  length  18  is  developed  from  this  noise  ..eries, 
and  is  used  to  predict  the  noise.  The  predicted  noise  series  are  shown 
as  the  third  graphs  from  the  top  In  each  part  of  Figure  4.  When  these 
predicted  noise  series  are  subtracted  from  the  input  noise  series  (top), 
one  obtains  the  filtered  times  series  (bottom).  By  comparing  the  bottom 
time  series  with  the  top  time  series,  one  can  determine  just  how  well  the 
noise  can  be  predicted.  Perfect  noise  prediction  would  manifest  Itself 
by  zeros  in  the  bottom  time  series.  Inspection  of  the  graphs  for  all  four 
channels  reveals  significant  noise  reduction  for  each  channel.  However, 
it  is  quite  obvious  that  the  noise  prediction  is  far  from  perfect. 

It  would  have  been  desirable  to  perform  further  experiments  with  the 
program.  Unfortunately,  because  of  the  untimely  termination  of  the  con> 
tract  this  was  not  possible.  Some  suggestions  for  further  work  is  con¬ 
tained  in  the  conclusions  and  recommendations  section. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  prediction-error  filter  technique  appears  to  offer  a  promising 
means  of  digitally  filtering  infrasonic  records  to  eliminate  some  of  the 
unwanted  noise  and  background.  In  this  technique,  a  sample  of  the  noise 
prior  to  an  Infrasonic  signal  is  used  to  develop  a  Wiener  prediction  fil¬ 
ter  for  predicting  the  noise  during  the  first  motion  (first  part  of  the 
event).  Application  of  the  filter  to  the  original  infrasonic  record  yields 
an  estimate  of  the  noise  before  and  during  the  first  motion.  When  this  is 
subtracted  from  the  original  record,  an  improved  estimate  of  the  first  motion 
is  obtained.  Results  based  upon  application  of  such  a  filter  to  several 
sets  of  infrasonic  records  -  including  two  samples  containing  infrasonic 
signals  and  one  sample  of  pure  noise  -  indicate  that  the  noise  does  have 
some  degree  of  predictability,  and  that  a  significant  amount  of  noise  re¬ 
duction  is  accomplished  by  the  filter.  However,  additional  experiments  are 
required  to  quantify  how  much  of  the  noise  Is  being  suppressed  and  to  deter¬ 
mine  optimum  values  for  the  parameters:  filter  length  and  prediction  span. 
These  experiments  can  best  be  performed  on  noise  samples,  since  one  can  com¬ 
pare  predicted  and  actual  values  of  the  I'oise.  A  sufficient  number  of 
experiments  should  be  performed  to  cover  the  various  types  of  noise  and 
background  that  can  occur  in  nature  and  to  permit  testing  of  various  com¬ 
binations  of  filter-length  and  prediction-span.  It  is  also  recommended  that 
an  averaging  procedure  be  added  to  the  technique.  This  would  permit  the 
averaging  of  the  filtered  time  series  from  all  the  channels  to  obtain  the 
best  single  estimate  of  the  first  motion. 

The  computer  program  MAXLKH  that  was  developed  under  the  contract  to 
perform  the  filtering  described  above  appears  to  be  operating  properly. 

Except  for  the  modifications  that  would  result  from  implementation  of  the 
recommended  experiments,  no  changes  are  presently  required  in  the  program. 

It  can  be  run  on  digitized  sets  of  infrasonic  data. 

Unfortunately,  because  of  the  untimely  termination  of  the  contract, 
we  were  not  able  to  pursue  the  other  aspect  of  our  research  program  -  the 
empirical  analysis  of  the  relationship  between  the  characteristics  of  a 
signal  and  the  meteorological  conditions  prevailing  along  the  path  of 
propagation  of  the  signal.  It  is  recommended  that  such  studies  be  under¬ 
taken  in  order  to  develop  a  technique  for  subtracting  out  the  effect  of  the 
atmosphere  on  the  propagation  of  an  infrasonic  signal.  Once  such  a  tech¬ 
nique  is  developed,  infrasonic  records  can  be  normalized  to  eliminate 
atmospheric  propagation  effects. 

Once  a  technique  for  eliminating  atmospheric  propagation  effects  is 
developed,  it  can  be  combined  with  a  noise  suppression  scheme  to  obtain 
an  optimum  estimate  of  the  actual  characteristics  of  the  original  infra¬ 
sonic  event.  Implementation  of  the  reconnendatlons  made  above  would  help 
to  achieve  this  goal. 
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